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Abstract 



We consider a pair of coupled non-linear partial differential equations de- 
scribing a biochemical model system. The Weiss-algorithm for the Painleve 
test, that has been succesfully used in mathemathical physics for the KdV- 
equation, Burgers equation, the sine-Gordon equation etc., is applied, and 
we find that the system possesses only the "conditional" Painleve property. 
We use the outcome of the analysis to construct an auto-Backlund transfor- 
mation, and we find a variety of one and two-parameter families of special 
solutions. 



1 Introduction 



Although the deeper reasons are still not completely understood, it is clear 
today that the Painleve property [1] of (systems of) ordinary or partial dif- 
ferential equations is closely related to the concept of integrability. The 
original approach of Ablowitz, Ramani and Segur [2], led to the conjecture 
that a non-linear partial differential equation (PDE) is integrable if every 
exact reduction to an ordinary differential equation (ODE) has the Painleve 
property, i.e. has no other movable singularities than poles. This version of 
the Painleve conjecture is well-suited to confirm already known integrability 
properties of a PDE, but is less helpfull when it comes to the act of finding 
new integrable systems, due to the problems of actually finding all the exact 
reductions. In this respect the later approach of Weiss, Tabor and Carnevale 
[3] seems to be much more usefull. In the simplest version of this approach 
a PDE is conjectured to be integrable if its solutions are singlevalued about 
movable " singularity- manifolds" : 

<f>(z 1 ,z 2 ,...,z n ) = 0, (1.1) 

where <f> is an arbitrary ("movable") analytic function. In other words a 
solution x(zi) to the PDE in question should have a Laurent-like expansion 
about the movable singular manifold (1.1): 

oo 

x(zt) = [^(zi^Y^M^M^Y , (i-2) 

3=0 

where p is a negative integer. If the number of arbitrary functions in the 
expansion (1.2) equals the order of the PDE, we have reason to believe that 
(1.2) is the generic expansion about = 0. If this is not the case we may 
have lost something and the system does not have the full Painleve property. 
Further details can be found in Refs. 1,3. 

This Weiss-approach to the Painleve property has been used to analyse 
a long list of the PDE's of mathematical physics like the KdV-equation, 
Burgers equation, sine-Gordon equation, modified KdV-equation, Boussinesq 
equation,... [3,4]. The method is of course also applicable to systems of 
coupled non-linear PDE's like the coupled KdV-equations (Hirota-Satsuma 
system) [5]: 

1 

x £ ~~ x I i^xXf Qyyy . 
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or the modified sine-Gordon system [5]: 

1 A x 

xt+~ x y+o - = o, 

2 2 y 

1 1 y 

2 2A x 

Systems of coupled non-linear PDE's also arise naturally in (mathemat- 
ical) biology and chemistry, for instance in models describing reactions and 
diffusion of different components in a medium. In this paper we will consider 
the following system of PDE's [6-10]: 

x t = x 2 y — Bx + Kx rr , 

y t = -x 2 y + Bx + Ky rr , (1.3) 

where B and K are constants. These equations are supposed to describe the 
evolution of concentrations x and y of 2 chemical components. Obviously K 
is a diffusion coefficient and the terms x 2 y — Bx represent the irreversible 
reactions. Such processes are certainly not expected to be integrable and 
we will indeed see that (1.3) does not have the full Painleve property (Sec. 
2). We note in passing that the system is closely related to the so-called 
Brusselator [6] describing the scheme of reactions: 

A^X, B + X + C, 
2X + Y -> 3X, X -> D, 

where A, B, C, D, X and Y are different kinds of molecules. 

The system of equations (1.3) has been investigated by several authors 
[6-10]. The steady state solutions (x t — yt — 0) were analysed in Ref. 7 
under the boundary condition x + y = const. In that case Eqs. (1.3) trivially 
separate and are completely solvable in terms of elliptic functions. In Refs. 
8 and 9 a system like (1.3), but with 2 different diffusion coefficients, was in- 
vestigated. It was assumed that one of the components diffuses very rapidly 
so that one of the diffusion coefficients could be taken to be infinite. The 
system was then reduced to a one-variable non-linear PDE and some trav- 
elling wave solutions were found. Finally we mention that a one-parameter 
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family of travelling wave solutions for the coupled system (1.3) was found in 
Ref. 8 under the boundary condition x + y = const. 

The paper is organized as follows: In section 2 we carry out the Painleve 
analysis of the system (1.3) using the Weiss- approach. We consider expan- 
sions of x and y about a movable singular manifold and find that the system 
possesses only the so-called conditional Painleve property. In section 3 we 
truncate the expansions of section 2 at the most singular term, which al- 
lows us to construct a two-parameter family of special solutions in terms of 
trigonometric functions. In section 4 we truncate the expansions of section 
2 at the "constant" term. This leads to an auto-Backlund transformation 
between 2 pairs of solutions. We then construct some more one and two- 
parameter families of special solutions, and in section 5 we finally give our 
conclusions and we outline how more general families of solutions can be 
"chased". 
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2 Painleve test 



In this section we perform the Painleve analysis of the system of coupled 
non-linear PDE's (1.3) with B and K as arbitrary positive constants (the 
positivity of B and K is of course not mandatory, at least not from a math- 
ematical point of view). 

The first step in the Weiss-algorithm [3] is to look for the dominant be- 
haviour about a movable singular manifold <p(r, t) = 0. Thus we write: 

x = U (r,t)[(f>(r,t)Y = U <l>P; p < 0, (2.1) 

y = V (r,t)l<P(r,t)r = V o r; *<0, (2.2) 

and balances the most singular terms after insertion in (1.3). This leads to 
the unique solution: 

p = a = -l, U = -V , (2.3) 

with: 

U = ±V2K<p r - (2.4) 

In the following we only consider the case of +- sign in (2.4). The other 
possibility corresponds to an overall change of sign of x and y. Clearly if 
(x, y) solves (1.3) also (— x, —y) is a solution. 

The next step is to look for the "resonances" [3], i.e. the orders in the 
expansions where arbitrary functions may appear. Keeping (2.3) and (2.4) 
in mind we write: 

x = ^2K(f) r <p- 1 + si0 ,_1 , (2.5) 
y = -V2K4> r (t)- 1 + s 2 0'-\ (2.6) 

where / is a non- negative integer and = Si(r,t);i = 1,2. After insertion in 
(1.3) and balancing of the most singular terms we get the matrix-equation: 

/ (I- 1)(Z — 2) -4 2 
^ 4 (Z-l)(Z-2)-2 

(2.7) 

with: 

det^= (I- 1)(Z — 2)[(Z — l)(Z-2) -6]. (2.8) 

For a given root of the polynomium (2.8) there are now as many arbitrary 
functions Sj as the multiplicity of that particular root. Solving the equation 
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detvA = we find the roots (—1,1,2,4). The root I = —1 corresponds to 
the arbitraryness of the location of the singular manifold (ft, while the roots 
I — 1, I — 2 and I = 4 are supposed to correspond to arbitrary functions at 
the orders 0°, (ft 1 and (ft 3 in the expansions of x and y about = 0. Note 
that the number of roots equals the order of the system (1.3) so until now 
everything works fine. 

The third and final step in the Weiss-algorithm [3] is to expand out to the 
highest resonance to make sure that no inconsistencies arise, i.e. we write: 



x 



= y/2K4> r <jT x + a 3 = a i( r ' *)> ( 2 - 9 ) 

3=0 

y = -V2KM- 1 + h = Pi(r,t). (2.10) 

3=0 

These expressions are inserted into (1.3) and we then balance the terms order 
by order in (/). It was demonstrated by Kruskal [11] that some simplifications 
arise in this process if one formally solves the equation 0(r, t) — for (say) 
r and then writes: 

<f>(r,t)=r + 4>(t), a j = a j (t), fa = /3j(t), (2.11) 

where ip is an arbitrary function. For our purposes of constructing explicit 
special solutions (sections 3,4) it is however necessary to use the general 
expansions (2.9) and (2.10). 

At the various orders we now get: 

V2K(2a -/%) = £- ^, (2.12) 

<Pr <Pr 

so here we get an arbitrary function (say) «o, as expected from the analysis 
of the resonances. 



a {2f3 - op) (f) rt B 
with another arbitrary function (say) a±. 



2K(2 ai - A) = — I ^ - 5 " X + K ^ ( 2 - 13 ) 

<Pr <Pr 9r <Pr 
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,0. 



a 2 = Fi(a , ai, 0; B, K) 
p 2 = F 2 (a ,a 1 ,(f); B,K) 



(2.14) 



where F\ and F 2 are certain complicated expressions in the arbitrary func- 
tions 0, «o, on and their derivatives, as well as in the "chemical" constants 
B and K. For convenience they are listed in the appendix. 



There is now an arbitrary function at this order if and only if the 2 right 
hand sides are equal. From the explicit expressions for G\ and G 2 given in 
the appendix it follows that this is not so (this is actually most easily seen by 
using the Kruskal Ansatz (2.11) and by keeping in mind the arbitraryness of 
(say) otQ and a\). It follows that (say) a% is arbitrary only if the singularity 
manifold 4> satisfies a certain constraint. Therefore, the system (1.3) does 
not have the full Painleve property but only the "conditional" one [12]. 

3 Truncation and special solutions 

In this and the following section we will extensively use the expressions among 
the alpha's and beta's obtained in section 2. In this section we look for special 
solutions to (1.3) obtained by truncation of the expansions (2.9) and (2.10) 
at the singular terms. Taking = $ = 0; % > we find: 



This can however only be a solution to (1.3) provided equations (2.12)-(2.15) 
are fulfilled. It is easily seen that quations (2.14) and (2.15) are trivially 
satisfied whereas (2.12) and (2.13) give the compatibility conditions: 




(2.15) 




(3.1) 




(3.2) 




rrr 



= 0. 



(3.3) 
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These 2 equations are integrated and consistently solved by: 

4>{r,t) = 



LU ' ci cos(W -^-r) + c 2 sin(W -^r) I f c 3 , 



2if 



where 01,02,03 are arbitrary constants. From (3.1) we then get: 



-%Bt 



e 2 



r) + L 1 sin( A /#r)) + 1 



(3.4) 



, . , . v / ^e-|^(L 1 cos(^r)-L 2 sin( v ^r)) 
x(r, i) = -y(r, t) = , —. ^ v - p= — T =- '-, (3.5) 



representing a two-parameter family (L x and L 2 being 2 arbitrary constants 
to be determined by the initial/boundary conditions) of solutions to the 
system (1.3). 



4 Auto-Backlund transformation and special 
solutions 

In this section we truncate the expansions (2.9) and (2.10) at order 0°, i.e. 
we take aij = fa — 0; % > 1, and look for solutions in the form: 

(:)^UH;:)- <«> 

As in section 3 this can of course only be a solution to (1.3) provided (2.12)- 
(2.15) are fulfilled. Equations (2.12) and (2.13) read: 

V2K(2a - /Wr = <k~ 3K<f>„, (4.2) 
a (2(3 -a ) = ^ + B-K^p. (4.3) 

<Pr <Pr 

After using (6.1) and (6.2) from the appendix, equation (2.14) leads to: 



a ot — ~ Bchq + Ka 



Orr j 



p ot = -a 2 (3 + Ba + Kf3 0rr , (4.4) 

while (2.15) is trivially fulfilled. (4.2)-(4.4) represents an overdetermined 
system of equations to be solved for a , /3 and 4>. This system is extremely 
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complicated in the general case but fortunately it is not so difficult to find 
special solutions. Note also that (4.4) has the same form as the original 
system (1.3), so if we can find solutions (a ,(3 ,(f)) to (4.2)-(4.4) it follows 
that (4.1) is an auto-Backlund transformation between the 2 pairs of solutions 
(x, y) and (a^, (3q) to equation (1.3). It means that if we find one pair of 
solutions we can in principle always generate new ones. 

Special solutions to the overdetermined system (4.2)-(4.4) can conve- 
niently be parametrized by the 2 functions [12]: 

C=£, V=^. (4.5) 

For simplicity we will now restrict ourselves by considering only constant 
C and V . It turns out that two different types of solutions are possible 
corresponding to V = and V ^ 0. In the case that V = and C = C is 
an arbitrary constant we find: 

<f>(r,t) = Cl (r + C t) + c 2 , (4.6) 

where Ci,c 2 are arbitrary constants. Equations (4.2)-(4.4) lead to: 

a (2/5 - a ) = B, (4.7) 
alfto - Ba = 0, 

that are solved by: 

a = p = iv^B, C = ±V2KB. (4.8) 
Using (4.1) we then get: 

\FlK 

x(r,t) = -y(r,t)±2y/B= 7= ± v^B, (4.9) 

representing a one-parameter family {L\ being an arbitrary constant) of so- 
lutions to (1.3). 

In the case that both V = V and C = C are arbitrary constants (Vo 7^ 0) 
we find instead: 

0(r, t) = Cl e Vo(r+Cot) + c 2 , (4.10) 
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where ci,C2 are arbitrary constants. In this case equations (4.2)-(4.4) lead 
to: 



V2K(2a - fa) = C - 3KV , 
a (2/9o - <*o) = V C + B- KV 2 , 



(4.11) 



«n/3o - -Ba = 0, 



Solving (4.11) for (a , (3 , C ) in terms of (B, K, V ) leads to the 3 possibilities: 

(4.12) 



(a ,A),C )= (-VwVo,^^,^- /vl ) . 



and: 



(a ,/?o,Co) = 



25 



v 2 -v^iO^iv/Wo 

where W = 2£? + -K^ 2 . From (4.1) we finally get: 



WqV2K ^ KV y/Wo \ 

-V2KV ±VW^ J ' 
(4.13) 



^/2KVoe Vo( - r+CotS) 
x(r, t) = -y(r, t) + a + ft = e y o(r+Cot) + ^ + «o, 



(4.14) 



with (ai ,/3 ,Co) given by one of (4.12), (4.13). This equation finally repre- 
sents a two-parameter family (Li and Vq being arbitrary constants) of solu- 
tions to (1.3). 



5 Conclusion 

In conclusion we have studied the pair of reaction-diffusion equations (1.3) 
using the Weiss-algorithm for the Painleve test. The system was found to 
possess only the conditional Painleve property. The results of the analysis 
however led to various kinds of special solutions obtained via truncations and 
auto-Backlund transformations. 

The special solutions we have constructed (3.5), (4.9) and (4.14) are all 
of the form x(r, i) + y(r, t) = const. Note that addition of the two equations 
in (3.1) leads to: 

(x + y) t = K(x + y)„, (5.1) 
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which is just the ordinary one- variable diffusion equation with the well-known 
general solution (see for instance reference 13): 



(x + y)(r, t) 



1 



— oo 



'OC 




(5.2) 



where: 



(x + y)(r,0) = F(r). 



(5.3) 



All our solutions therefore correspond to the boundary condition F(r) = 
const., and could in principle have been obtained by applying the Weiss- 
algorithm to a one-variable non-linear react ion- diffusion equation obtained 
from (3.1) by separating x and y = —x + const, from the beginning. 

The advantage of our more general approach where (in the notation of 
(5.2)) F(r) is not fixed from the beginning , is of course that we can now 
use the results of sections 2-4 to look for more general families of solutions 
with non-constant F(r): We can either return to equation (4.5) and look for 
solutions to (4.2)-(4.4) with non-constant V and/or C, or we can continue 
along the road of sections 2 and 3 and look at expansions truncated at higher 
orders in (j). That is however out of the scope of this paper. 

Acknowledgements: I would like to thank J.B.Pedersen for a discussion 
on general aspects of diffusion and reactions. 
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6 Appendix 



In this appendix we list the functions (Fi,F 2 ) appearing in (2.14) and the 
functions (Gi,^) appearing in (2.15). 
The functions F\ and F 2 are given by: 

4K&F! = p ot - (3 l( f> t + a 2 (3 - Ba - K/3 0r r - 2KM r 

- K(3i<j) rr + 2V2K(p r (a f3i + ai(3 - a ai), (6.1) 



4K$F 2 = (2a + Po) t + (2«i + A)0t - K(2a + /3 ) rr 

- 2K(2 ai + ft) r r + Ba - a 2 (3 - K(2a x + 

- 2KV2K(p r (aof3 1 + a^fa - a ai). (6.2) 

Using (2.12) and (2.13) it is straightforward to express F\ and F 2 in terms of 
(say) the arbitrary functions (a , eti, <fi) and their derivatives as well as the 
"chemical" constants (B,K). 

The functions G\ and G 2 are given by: 

AK^d = 2a lt + 4a 2 <f) t + 2V2K4> r (2a a 2 + aj) + 2Bu x 

- Ap (V2K(j) r a 2 + a ai) - 2/3 1 (2v / 2i^> r ai + a 2 ) 

— 4\/2K/3 2 (p r ao — 2K[a\ rr + 4a 2r (p r + 2a 2 (j) rr ), (6.3) 



AK<j) 2 .G 2 = (3 lt + 2p 2 4> t - V2K(P r (2a a 2 + a 2 ) - Bon 

+ 2(3 (V2K(j) r a 2 + a «i) + A(2v / 2K0 r ai + a 2 ) 

+ 2V2K(3 2 (f) r ao - K{(3 lrr + 4p 2r (f> r + 2/3 2 rr ). (6.4) 

Using (6.1), (6.2) and (2.12), (2.13) we can again express the right hand sides 
in terms of (say) the arbitrary functions (cto,ai,4>) and their derivatives as 
well as the "chemical" constants (B, K). It can then be verified that G\ ^ G 2 , 
so that both a 3 and (3% are fixed by (2.15), i.e. there is no arbitrary function 
at this order. 
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